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.  PROPAGATION  OF  SOUND  IN  THE  HOMOGENOUS  LAYERED  MODEL 


The  transmission  of  sound  in  the  ocean  is  dependent  on 
several  environmental  factors.  The  most  important  of  these 
include  the  depth,  the  nature  of  the  bottom,  sound  velocity 
structure  within  the  ocean,  and  the  shape  of  ocean  bottom 
and  surface.  The  purpose  of  this  experiment  was  to  investi¬ 
gate  the  propagation  of  sound  in  a  two  layered  laboratory 
model,  composed  of  sand  and  water.  In  this  study  the  bottom 
and  surface  of  the  model  were  made  flat  and  experimental 
results  were  compared  with  a  simple  theoretical  model. 


II.  PROPAGATION  OF  SOUND  IN  SHALLOW  WATER 


A.  INTRODUCTION 

The  propagation  of  sound  in  shallow  water  is  a  compli¬ 
cated  problem,  involving  multiple  reflections  of  the  sound 
between  the  surface  and  the  bottom,  so  that  there  is  consid¬ 
erable  phase  interference  among  the  various  multipaths.  At 
relatively  long  ranges,  the  sound  field  is  mainly  due  to  the 
multiple  reflections.  The  transmission  of  sound  in  the 
shallow  water  also  depends  on  several  environmental  factors: 
nature  and  slope  of  the  bottom,  depth  of  water  relative  to 
wave  length,  roughness  of  surface  and  bottom  and  temperature 
gradient.  As  a  consequence  of  the  multiplicity  of  reflec¬ 
tions  from  the  surface  and  the  bottom  and  the  environmental 
factors,  the  predictions  of  sound  transmission  in  shallow 
water  can  be  much  more  complicated  than  for  deep  water.  Ray 
theory  can  be  used  for  the  sound  field  at  short  range  when 
the  depth  of  the  water  is  relatively  large.  But,  ray  theory 
has  limited  applications.  It  is  inapplicable  in  a  shadow 
zone.  Ray  theory  also  cannot  be  used  for  low  frequencies 
when  the  wavelength  of  sound  becomes  comparable  with  the 
vertical  scale  of  sound  velocity  variation.  Therefore,  in 
practice  one  ultimately  has  to  resort  to  a  wave  solution  of 
the  problem. 

B.  NORMAL  MODES  IN  TWO  Hv.  '  '’OUS  LAYERS 

Normal  mode  theory  provxu  an  exact  solution  to  the 
wave  equation  and  is  also  use  .1  for  dealing  with  depth 
dependent  sound  profile  partic.  Larly  those  leading  to 
the  trapping  of  sound.  There  a.=  various  approaches  to 
the  normal  mode  theory  but  we  will  use  the 


tan  Xn  -  —  X„ 

bHVi^(5j 


(2.55) 


tan  Xn 
X„ 


-  1 _ 

b  h/  _  Xyi  (2.56) 

c.*“  c/  H2 


When  C2  >  C,  which  is  the  case  under  consideration, 
equation  (2.56)  possesses  real  roots  for  k„  and  X^,  ,  the 
later  varying  in  the  range  (n  -  1/2)  IT  £  X^£  n  <Y  .  On  the 
other  hand,  when  C2  <  C,  ,  the  roots  k n  are  complex  numbers 
with  negative  imaginary  parts.  In  the  later  case,  there¬ 
fore,  the  factor  exp(-jk^r)  in  equation  (2.52)  implies 
horizontal  attenuation.  In  the  case  of  a  fast  bottom,  there 
exist  solutions  which  suffer  no  horizontal  damping.  This 
condition  stems,  of  course,  from  the  fact  that  for  angles  of 
incidence  (measured  from  the  normal)  greater  than  critical 
angle,  no  power  is  transmitted  into  the  bottom.  Figures 
2.2,  2.3,  and  2.4  shows  the  z  dependence  of  the  first  three 
normal  modes  for  50  kHz,  120  kHz,  160  kHz  respectively. 
Figure  2.5  shows  the  r  dependence  of  the  sound  pressure. 
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Solution  to  this  homogeneous  equation  is 


R„(r>  =  -j7TZ„(d)H0(k„r) 


(2.51) 


where  for  knr  >>  1 


(2) , 

H  (k 


cil 

The  Hankel  function  H0  ( k  ^  r)  describes  the  solution  in 
the  r  direction. 

Therefore,  the  solution  for  the  sound  trapped  in  layer  1 
(0  <  Z  <  H)  is 


P<r,z,t) 


A„sin(kd)sin(k  z 

n  zn  zn 


)e/r«'-V  + 

(2.52) 


Equation  (2.52) 
of  Z  and  d  and  is 
The  eigenvalues  k2„ 
equation  (2.25) 


is  invariant  under  the  interchange 
valid  within  the  region  0  <  Z  <  H  . 
can  be  obtained  by  eliminating  Bn  from 


tan  k  H 

zn 


.Aid 
A  B, 


(2.53) 


Let 


X 


n 


k  H 

zn 
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j2 

A 


(2.54) 


The  depth  eigenfunctions  now  form  a  complete  set  of 
orthonormal  functions. 

Mow  we  have  to  evaluate  the  Green's  function  for  the 
inhomogeneous  wave  equation 


S(r)S(z-d) 

(2.46) 


2  Sir)  S(z-d) 
r 


(2.47) 


From  equation  (2.13),  the  second  term  in  the  square 
bracket  is  zero;  therefore,  equation  (2.47)  can  be  rewritten 
as 


Izrldfrd&A  k  R 
n  "L  r  d  rV  d  r  >  n~n 


_  2_  Sir )  6(2  -d  ) 

r  (2.48) 


Multiply  both  sides  of  equation  (2.43)  by  Zwand 
integrate  from  0  to  co  gives 


y  f  1  iL(  r  dJir)  +  k  R  1  r  Z„Z„d  z  -  _2_&(r)f&(  z-d)Z„dz 
V1  r  d  rv  (jrJ  n~n\J  "  r  7 


=  n 


iX('rdfi»'\+  k*R  = 

r  dr'  dr  /  + 


2  6(r)Z  (d) 
r  * 


For  m 


Integrating  and  substituting  the  limit  we  get 


and 


or 


I 


H  2 

Z  dz 


in 


=  ±  f  H  _  J_  sin  k  ■  H  cos  k  -  H  1 

W  l  T  m  in  j  ( 2 . 42  ) 


'in 


r03  2  s 

J  ^rfi1  =  f 

*L/  X  i 


(  sin2  k  H-e^  )e*^dz 


in 


H 


H 


f 


co  2 

Z.  dz 


in 


1  Sin2k  H 

—  in 


H 


2& 


(2.43) 


(2.44) 


Substituting  equation  (2.42)  and  equation  (2.44)  into 
equation  (2.40)  gives 


2  k 


‘in 


k  H  - coS(k  H)sin(k nH)  -/fifsinfk  H) -tan(k nH) 

in  m  m  j  in  in 


in  in 

(2. AS) 


This  is  only  possible  if 


i  .  7  A.  =  o 
A 


or 


7 


1 

(2.36) 


Therefore,  equation  (2.29)  becomes 


OS  2  2  H 

W<  Z<  oo  (2.37) 

We  will  normalize  the  depth  eigenfunctions  in  order  to 
evaluate  An  ,  consider  the  m  =  n  case 


V 

Zn{ 2 


)d  z 


(2.38) 


(2.39) 


(2.40) 
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interchange  the  subscript  m  and  n  in  equation  (2.30)  and 
take  the  difference  of  the  two  equations  to  get 


2  .2 


cL(Z  dZp  __  Z  =  (  k„  -  k  ') Z  z 
dz'  dz  )  \  n  m)  m  n 


(2.31) 


Integrate  both  sides  of  equation  (2.31)  from  0  to 


2  2  .  r* 


(2.32) 


The  orthogonality  condition  requires  that,  for  m  ^  n, 
the  R.H.S  of  the  equation  (2.32)  must  be  zero.  Therefore 


H  2 

(  dl.n_2,„dZ1«\dz+7  [d_fZ  dZ„  z  dZ„\dz  -  0 

Jdz^’-dz  dT  '  +  J  dz(  "d~“  "dT^ 


(2.33) 


m  _  2 


[  Z,„  dJL.*  _  Z,„  dZ_,«j  4.  T  [  dZ in  _  Zln  =  o 

L  dZ  HZ  l  I  H  7  d  z  JH 


(2.34) 


The  application  of  limits  as  well  as  boundary  conditions 


gives 


[ 1  -7  A  ] [  ZjH>(dji  -  z»tH,(ffr’)  ]  =  0 


(2.35) 


I 


which  gives 


B„ 


(2.26) 


The  depth  function,  therefore,  can  now  be  written  as 


Z,  Cz)  =  A  sin  k • z 

in  n  zn 


0<2  <  /f  (2.27) 


z 


A  sin  k-  H ■ 

n  zn 


ep'(7'") 


H  c  2 


00  (2.28) 


These  depth  functions  have  to  be  orthogonal.  Since  there 
is  a  discontinuity  in  the  slope  at  the  boundary  Z  =  H,  we 
need  separate  orthogonality  constant  for  each  layer.  We 
assume  that  the  constant  is  unity  for  the  layer  1  and  7S  for 
the  layer  2. 


z„(z> 


Z,  (z)  o<z<« 

1/) 

y  Z2„(z) 


In  order  to  evaluate  75  multiply  equation  (2.13)  by  Zw 


2 

Z  dZ 
"dl? 


+ 


2  2 
(k  -  k 


)Z  Z 

n  n  m 


=  0 


(2.30) 
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the  boundary  conditions  at  the  surface  at  Z  = 
be  sinusoidal  function  of  depth.  Therefore 


His  to 


z 


in 


A  Si n  k  2 

n  zn 


0-5.Z-S  H 


z 

2D 


H<Z-<oO 


ft 


Waveguide  propagation  requires  that  kM  and  must 
hence 


"2k 

C< 


Co 


Application  of  boundary  Conditions  at  Z  =  H  results 


Ansin  k •  h 

zn 


A„  J_  k  cos  k  h 
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_  j_  BBe 

n  &  n 


(2.20) 


(2.21) 


(2.22) 

be  real 


(2.23) 

in 


(2.24) 
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(2.25) 


/?VdZ  JH 


yO  v  d  Z 


.  lb) 
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l'^Z2fl(z)  =  0 

2—0 o 


Rewrite  equation  (2.13)  in  the  form 


2  .2 


^.+(K  -M2- 


=  0 


or 


a-7  2 

d  Z  n  ,  j<  2 


=  0 


where 


In  order  to  satisfy  the  boundary  condition 
goes  to  zero  when  Z  — ►  oo  ,  the  solution  Zt„ 
exponentially  with  depth.  The  only  way  that  Zm 


(2.17) 


(2.18) 


(2.19) 

the  pressure 
must  decay 
can  satisfy 


Substitute  equation  (2.9)  into  equation  (2.8)  to  get 


f  Y*  Z  J_  (  r  dBjj)  ,  Y*  8  sL?*,  +  Z  B  Z  k  1  = 
L  „  %r  3 r ^  ar  »  fl0P+  "  -1 
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(2.10) 


or 


r  j_jl 

d_( 

'  r  dRpV 

1  +[ 

L  r 

dr' 

•  dr  J 

=  0 


(2.11) 


This  is  only  possible  if  each  part  in  the  square 
brackets  in  equation  (2.11)  is  equal  to  a  constant  hence  it 
can  be  written  as 


1±  d_/rd8*A  =  _  kf 

R  r  drV  dr  J 


(2.12) 


X  j-  Z,  +  k 

z„di^ 


=  K 


(2.13) 


where  kn  is  separation  constant. 

The  boundary  conditions  on  the  depth  functions  are 


where  k  is  wave  number  defined  by 


k  = 


CO 

c 


(2.5) 


The  delta  function  for  the  point  source  at  r  =  d  is 


S(r-x,) 


_  1  S(  r)  6(Z  -d  ) 
r 


(2.6) 


In  cylindrical  co-ordinates  equation  (2.4)  can  be  written  as 


fi£_(ra_\  af  k2  ]P  =  _  2_  6(r)c5(  z-d)  ey"' 

L  r  aA  ar/  3z2  +  J  r 


(2.7) 


This  is  an  inhomogeneous  equation.  First,  let  us 
consider  the  homogeneous  equation. 


r  LJL( 

r3  \ 

L  r  an 

dr  ' 

+ 


3  +  k  TP  =  0 


(2.8) 


and  use  the  separation  of  variables  technique. 

We  assume  that  the  sound  energy  is  trapped  in  the  upper 
layer  and  solution  to  equation  (2.8)  can  be  represented  by 
a  set  eigenfunctions. 


P(r,z,t) 


Z  R(r)Z(z) 

n 


n 


(2.9) 


14 


f 


4 


> 


I 

t 

k' 

S 


If 


V  p  -  1  a2p  =  -  6 

-  PaF 

where  <T  is  a  source  function  and  is  defined  by 

or  =  6Q 

at 


> 


i 


Where  Q  is  source  strength 

In  our  solution,  we  assume  that  a  point  source  radiating 
a  single  frequency  can  be  described  mathematically  by 


Cj 


4TTA  S(  r -4) 


(2.2) 


Where  S(r  -  r.  )  is  the  Dirac  delta  function,  r  is  the 
position  vector.  A  is  the  amplitude  factor  which  is 
assumed  to  be  unity,  and  co  is  the  angular  frequency  of 
the  source.  Solution  to  the  wave  equation  is  of  the  form 


P(r,t)  =  P(r)e/W 

~  ~  (2.3) 


Substituting  equation  (2.2)  and  equation  (2.3)  into 
equation  (2.1)  we  get  the  helmholtz  equation 


Vp 


2 

k  P 


=  -  41T Sir-  r  ) 


(2.4) 


» 


1 
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eigenfunction-expansion  technique,  which  is  relatively 
simple.  In  this  presentation  we  will  follow  mainly  [Ref.  1] 
and  [ Ref .  2 ) . 

Consider  the  case  of  two  homogeneous  layers  with  a 
pressure  release  surface  above  as  shown  in  Figure  2.1  . 


Figure  2.1  Two  homogeneous  layers  with  pressure  release 
surface . 

Let  and  f ,  be  the  constant  densities  and  sound 
speeds  in  layer  1  and  2  respectively.  Assume  there  is  a 
point  source  at  z  =  d  in  layer  1.  The  wave  equation  for  a 
homogeneous  medium  can  be  written  as 


Figure  2.4  Pressure  amplitude  vs  depth  for  third  mode 


tude  as  a  function 


C.  CUT  -  OFF  FREQUENCY 


Let  us  determine  the  cut-off  frequencies  of  the  normal 
modes.  It  follows  from  equation  (2.56)  that 


It  is  clear  from  the  above  equation  that  the  nth  mode 
cannot  support  propagation  until 


(2.58) 


Therefore,  the  limiting  frequency  below  which  transmis¬ 
sion  of  the  nth  mode  cannot  take  place  is  given  by 


f  -  (2n  - 1  )c, 

n  ~ 

where  n  =  1,2,3, - 

This  is  an  important  relation  between  the  layer  thick¬ 
ness  H  and  the  cut-off  frequency  f  .  It  is  shown  in  Figure 
2.6  for  the  first  mode. 


4 HJ  1  -  c? 

Cl 


(2.59) 


Cut-off  frequency  as  function  of  layer  thickness 


D.  PHASE  SPEED 

The  factor  exp  j  (tot-k^r*  1^4)  in  equation  (2.52)  allows 
an  obvious  interpretation  of  kM,  namely  kn  =  Cpn  where 
CPn  is  the  phase  speed  of  the  nth  mode.  The  phase  speed  is 
the  speed  that  the  wave  front  (surface  of  constant  phase) 
propagates  in  r  direction.  The  phase  speed  at  the  cut-off 
frequency  has  the  value  of  the  speed  of  sound  in  the 
bottom  and  decreases  continuously  with  the  increasing 
frequency  towards  the  value  of  sound  speed  in  the  upper 
layer.  From  equation  (2.19)  we  can  write  (Ref.  3] 


2  2  2 
K  =  tei  -  K. 


Multiply  both  sides  by  H  gives 


(2.60) 


2  2 


,  j.  2  2  2 

=r  (  wH  \  _  k  H 

VC,/  zn 


(2.61) 


cj  H 


2  2  2  2  2 
^(kH+k^H  ) 


c/h2  = 


2  2  2 

C,  (  \  } 


(2.62) 


where 


2 


Therefore  equation  (2.57)  can  be  written  as 


ton  Xn 
X„ 


or 


tan  Xn 
X 


where 


6  =  Al± 


From  equation  (2.65) 


CX 

n 


X^V-T+  5  scot*(  X  ) 
6  n 


(2.63) 


(2.64) 


(2-65) 


(2.66) 


(2.67) 


(2.68) 
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We  know  that  C(n=  knH  where  kn  -  to/Cvr,  . 
Cp„  is  the  phase  speed  of  the  nth  mode. 


Therefore 


C  = 

*r> 


(2.69) 


From  equation  (2.62)  we  have 


i  I  r 

*  +  X 
n  n 


2.70) 


Hence 


=  C,V1  +('j<n 


(XJL) 

\  <X„  ) 


(2.71) 


Substituting  the  value  of  Xn/0<n  from  equation  (2.68)  we  get 


C  = 

PK» 


<V  1  + 


1  +  6  cot(Xn5 


(2.72) 


( n- 1/2 )  IT  <  X„  <  n  IT  for  the  nth  mode,  and  in  this  case 
k„  and  phase  speed  are  real  and  propagation  of  the  modes 
proceed  without  damping.  X„  is  a  function  of  angular 
frequency  and  the  frequency  corresponding  to  any  given 


value  of  Xn  can  be  obtained. 
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7  E 


cp  H 
27TC, 


2  7T6 


=  Cp. 

2  7TC, 

r  *  7*  * 

1  +  €  +  6-cot  X„ 


(2.73) 


When  Xn  is  near  its  upper  limit  of  n'iT'  ,then  cot  X  -*•  oO. 
Therefore  ft  is  very  large  and  from  equation  (2.72) 


Hence  in  the  limit  of  very  high  frequency,  phase  speed 
approaches  the  sound  speed  CH .  When  Xn  is  near  its  lower 
limits  (n-l/2)'Tf  then  cot  X  — *  0  . 

and  Cfn=  C,[  1  +  6-j1  =  c2 

The  phase  speed,  therefore,  varies  from  the  value  Cx  at 
the  critical  frequency  down  to  C,  in  the  limit  of  very  high 
frequency.  Figure  (2.7)  shows  the  phase  speed  as  a  function 
of  frequency  for  first  three  modes. 


E.  GROUP  SPEED 

In  a  medium  in  which  the  phase  speed  varies  with 
frequency  and  if  the  pressure  pulse  from  a  point  source  is 
not  a  single  frequency  but  covers  a  moderately  broad 
spectrum,  then  the  presence  of  the  particular  frequency 
component  of  the  spectrum  at  a  large  range  would  occur  at  a 
time  associated  with  the  group  speed  (speed  with  which  the 
energy  associated  with  a  particular  frequency  will  travel). 
At  frequencies  co  for  which  the  phase  is  stationary, 
mutual  interference  will  be  negligible  and  these 
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frequencies  will  be  dominant  at  the  prescribed  values 
of  t  and  r.  The  values  of  t  and  r  for  which  the  phase 
is  stationary  for  a  given  frequency,  therefore,  determine 
the  rate  of  propagation  of  this  frequency  in  the  mutually 
interfering  train  of  sinusoidal  waves,  [Ref.  3],  called 
group  speed  . 


but 


Therefore 


f(w,r,t)  - 


cot  -  k  («)  r  4-  7T 
n  4 


(2.74) 


d_L  =  t  -  r  d_kn_  =  0 

dcu  d  co 


(2.75) 


c  „  =  r 

1n  T 


d  co 
d  k„ 


(2.76) 


CO 


= 

ff  Pn 


d(knC,„)  = 

d  k„ 


k  d  cPr1 
nd  k„ 


(2.77) 


From  equation  (2.63)  we  have  k„  =  (X„/ H  .  Substituting  this 
value  in  equation  (2.77),  we  get 
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(2.78) 


C 


CPn  + 

dc<n 


From  equation  (2.71)  we  can  write 


From  equation  (2.68)  we  get 


d  X/j  _  6  /  1  +  S2cot2  x'n 

d«n  0(X„) 


(2.80) 


And  then  substitution  of  equation  (2.80)  and  equation  (2.79) 
into  equation  (2.78)  yields 


where 


(2.81) 


1  +  b cot‘?Xn  - 


8  xn  cos  Xn// s\n\  J 


(2.82) 
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Figure  (2.7)  shows  the  group  speed  as  a  function  of 
frequency  for  first  three  modes.  It  may  be  noted  that  group 
speed  is  always  less  than  the  phase  speed.  The  general 
features  for  each  mode  are,  [Ref.  2]. 

1. The  first  arrival  would  be  nearly  sinusoidal  waves 

of  frequency  f  n  ,  where  fM  denotes  the  limiting 

frequency  for  the  n  th  mode.  The  wave  would  arrive 
at  t  =  r/C2  travelling  with  C9yi  =  C2  .  This  portion 
of  the  signal  is  called  ground  wave.  As  the  time 
progresses,  the  frequency  and  amplitude  in  the  ground 
wave  would  increase. 

2.  At  time  t  =  r/C1  a  new  train  of  high  frequency 

waves,  due  to  the  right  hand  branch  of  the  group 
velocity  curve,  would  suddenly  be  superimposed  on 
the  ground  wave  travelling  with  speed  =  C,  . 

This  new  high  frequency  wave  is  designated  as  water 
wave.  The  frequency  in  the  water  wave  would 
decrease  as  time  progresses  while  in  the  ground 
wave  the  frequency  would  continue  to  increase. 

3. Still  later  fa  and  fw  would  approach  each  other 
at  t=r/Cg0  ,  where  Cg0  denote  the  minimum  group  speed 
they  would  coincide.  The  pressure  then  would 
consists  of  a  single  frequency  fA  and  this  portion 
of  the  signal  is  called  Airy  phase.  The  sequence  of 
these  events  is  illustrated  in  Figure  2.8  . 
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TABLE  V 


Sound 

f 

pressure  as 
-  35  kHz,  n 

function 
=  1,  L  = 

of  depth: 

1.5  m 

' 

No 

Depth ( cm ) 

Volts 

Normalized 

1 

0.0 

0.012 

0. 12 

2 

0.5 

0 . 028 

0.28 

3 

1.0 

0 . 056 

0.56 

4 

1.5 

0.081 

0.81 

5 

2.0 

0.097 

0 . 97 

6 

2 . 5 

0.1 

1 . 0 

7 

3 . 0 

0.097 

0.97 

8 

3 . 5 

0.083 

0.83 

9 

4.0 

0 . 067 

0 . 67 

10 

4.5 

0.052 

0.52 

... 

11 

5.0 

0.042 

0.42 

12 

5.5 

0.041 

0.41 

13 

6.0 

0.038 

0.38 

14 

6.5 

0.034 

0.34 

15 

7.0 

0.028 

0.28 

16 

7 . 5 

0.023 

0.23 

17 

8.0 

0.018 

0. 18 

18 

8.5 

0.015 

0. 15 

* 

19 

9.0 

0.016 

0 . 16 

53 


E. 


VALIDATION 


To  insure  that  the  sound  fields  set  up  in  the  tank  would 
be  good  approximations  of  the  sound  fields  in  the  absence  of 
sides  and  bottom  (e.g,  a  shallow  water  environment  of  infi¬ 
nite  horizontal  extent),  measurement  were  made  on  the  sound 
field  set  up  in  the  water  and  in  the  sand  for  a  source  in 
the  water  layer. 

1 .  Sound  Pressure  as  a  Function  of  Receiver  Depth 

Figure  3.5  shows  the  circuit  arrangement  for  the 
measurement  of  the  variation  of  pressure  amplitude  with 
range  and  depth.  Continuous  waves  of  the  desired  frequency 
were  generated  by  a  frequency  generator  and  were  applied  to 
the  source  through  an  amplifier.  The  source  was  placed  at 
mid  depth  in  the  water  layer  which  was  2.7  cm  thick.  The 
signal  was  received  by  the  receiver  and  was  passed  through 
the  bandpass  filter  to  eliminate  noise.  The  output  from  the 
bandpass  filter  was  amplified  and  displayed  on  the  oscillo¬ 
scope.  To  determine  the  variations  of  pressure  amplitude  as 
a  function  of  receiver  depth,  the  transmitter  was  kept  at  a 
constant  depth  in  the  water  layer  whr  le  the  measurements 
were  made  with  the  receiver  at  a  fixed  range  beginning  from 
the  top  surface  of  the  water  layer  down  into  the  bottom 
layer  in  steps  of  0.5  centimeter.  Three  frequencies  between 
the  cut-offs  of  two  first  modes  were  selected  for  measure¬ 
ments.  The  data  are  tabulated  in  Tables  V,  VI,  and  VII  and 
plotted  in  Figures  3.6,  3.7,  and  3.8  .  The  experimental  and 
theoretical  results  agree  satisfactorily  within  the  water 
layer.  Near  the  cut-off  frequency  of  the  propagating  mode, 
the  experimental  amplitude  appears  to  decay  more  rapidly 
with  depth  in  the  bottom  layer  than  theory  predict.  Decay  is 
slower  for  intermediate  and  high  frequencies. 


TABLE  IV 

Voltage  response  and 

sensitivity 

- 

Fre . 

Sensitivity 

Voltage  Response 

(kHz) 

dB ( re  V/Pa) 

dB ( re  Pa/V  §  1  m) 

20 

-76.37 

5.91 

30 

-68.34 

0.38. 

40 

-70.96 

1.06 

50 

-63.43 

6.84 

60 

-65.78 

8.83 

70 

-59.36 

12 . 50 

80 

-59.91 

11.03 

90 

-53.95 

15.39 

100 

-54.24 

16.18 

110 

-48.96 

20.32 

120 

-61.86 

26.87 

130 

-63 . 33 

27.94 

140 

-68.10 

29 . 69 

..... 

150 

-63 . 68 

29.07 

•  V 

TABLE  III 


Experimental  measurements  for 
calibrations  of  transducers 


Ere . 
(kHz) 

Vr 

(volt) 

Vx 

(volt) 

U 

(amp) 

20 

0.00046 

0.0015 

0.0051 

30 

0.00073 

0.002 

0.0075 

40 

0.0012 

0.0016 

0.01 

50 

0.00175 

0.0074 

0.0127 

60 

0.00255 

0.0071 

0.0155 

70 

0.003 

0 . 0227 

0.0185 

80 

0.0046 

0.018 

0.0215 

90 

0.006 

0.059 

0.025 

100 

0.0075 

0.0625 

0 . 0275 

110 

0.01 

0.185 

0.0305 

120 

0.0123 

0.089 

0.0315 

130 

0.016 

0.085 

0.034 

140 

0.015 

0.06 

0.037 

150 

0.0145  ' 

0.093 

0.043 

V* 

( volt ) 
0.00036 
0.006 
0.012 
0.034 
0.044 
0 . 099 
0.224 
0.46 
0.62 
1 .  15 
0.17 
0.193 
0.1 
0.21 


Distance  between  source  and  receiver 
Applied  voltage 


1  m . 

5  volts 


Figure  3. 


k -  r  =  1  m. 


1 

=P2I 

2 

/.  =x«\ 

3 

7  =  Transmitter  (  Source ) 
R  =  LC-10 

X  3  Receiver 


2  Calibration  of  transducers  by  reciprocity  methc 
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D.  VOLTAGE  RESPONSE  AND  SENSITIVITY  OF  TRANSDUCERS 


The  basic  characteristics  of  the  source  and  receiver 
were  not  known,  therefore  the  voltage  response  of  the  trans¬ 
mitter  and  the  sensitivity  of  the  receiver  were  measured  by 
the  reciprocity  method  (Figure  3.2).  The  left  hand  side  of 
the  Figure  indicates  the  transducers  being  used  as  transmit¬ 
ters  and  the  right  hand  side  as  receivers.  Experimentally 
measured  data  are  shown  in  Table  III  and  calculated  values 
of  voltage  response  and  sensitivity  are  shown  in  Table  IV. 
The  voltage  response  of  the  transmitter  is  plotted  in  Figure 
3.3  and  sensitivity  of  the  receiver  in  Figure  3.4  . 


TABLE  II 


O 

Speed  of  sound  in  water  (Temperature  =17  C) 


No. 

Time 

Division 

Scale 
( s/div) 

50 

Distance 

(cm) 

C 

(m/s ) 

1 

3.8 

27 . 5 

1447 . 3 

2 

4.5 

50 

33.0 

1466.67 

3 

6.0 

50 

43 . 5 

1450.0 

4 

7.05 

50 

51.5 

1461 . 0 

5 

8.8 

50 

64.5 

1465 .91 

6 

5.0 

100 

72.2 

1444.0 

7 

5.7 

100 

83 . 7 

1468.42 

8 

6.4 

100 

93 . 7 

1464.06 

Average 

Speed  of  Sound  = 

1458.43  m/s 

Standard 

Deviation 

= 

9.1 

Speed  of 

Sound  from 

Theory  = 

1472.39  m/s 

Percent  Deviation 


1472.39  m/s 
0.95  % 


TABLE  I 

Speed  of  sound  in  sand 


No. 

Time 

Division 

Scale 
(  A  s/div) 

Di stance 
(cm) 

C 

(m/s ) 

1 

8.1 

50 

67.8 

1674.07 

2 

6.4 

50 

54.3 

1696 . 88 

3 

5.5 

50 

45.3 

1647 .27 

4 

4.6 

50 

38.0 

1652 . 17 

5 

6.3 

20 

21.1 

1674. 6 

6 

8.6 

10 

14.2 

1651 .  16 

Average  Speed 

= 

1666.03 

m/s 

Standard  Deviation 

= 

17 . 58 

Percent  Deviation 

= 

1.06 

% 

Measurement  For  Density  of  Sand 
Mass  of  the  Container  =  542  gm 

Mass  of  1000  cm3  of  sand  +  Container  =  2650.3  gm 
Density  =  2108.3  kg/m3 
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Figur  3.1  Measurement  of  speed  of  sound  using  1  MHz  tone  burst. 


deep  into  the  sand  with  a  centrifugal  pump  at  moderate  pres¬ 
sure.  This  process  was  repeated  several  times.  The  top 
surface  of  the  sand  was  leveled  with  reference  to  the  water 
surface.  This  leveled  surface  of  sand  was  horizonally  flat 
within  +  1  mm.  The  layer  of  water  was  very  shallow  and  the 
thickness  of  this  layer  changes  with  time  due  to  evaporation 
of  the  water.  To  avoid  this,  the  tank  was  equiped  with  an 
automatic  water  level  control  which  maintains  the  desired 
level  of  water. 

2.  The  frequencies  used  in  this  work  ranged  from  25  kHz 
to  100  kHz.  The  source  transducer  was  built  from  a  rectan¬ 
gular  crystal  of  piezoelectric  material  2.98  cm  long  and 
1.12  cm  wide.  The  directivity  index  of  the  transducer  at  50 
kHz  is  5.75  dB.  The  receiver  was  made  from  a  4  mm  by  4  mm 
piezoceramic  cylinder. 

C.  MEASUREMENT  OF  SPEED  OF  SOUND  AND  DENSITY 

Figure  3.1  shows  the  equipment  and  the  circuit  used  to 
measure  the  speed  of  sound  in  water- saturated  sand  and  in 
water.  Two  LC-5  transducers  were  used  as  source  and 
receiver.  A  high  frequency  tone  burst  of  the  order  of  one 
MHz  was  used  and  the  speed  of  sound  was  calculated  by 
direct  measurement  of  time-of-f light  from  the  oscilloscope 
and  the  distance  between  the  source  and  the  receiver. 
Experimentally  determined  values  for  the  speed  of  sound  in 
the  sand  and  in  water  are  shown  in  Table  I  and  Table  II.  The 
standard  deviation  of  the  measurements  is  an  indication  of 
the  precision  of  the  measurement  (approximately  1  percent). 
The  difference  between  the  measured  speed  of  sound  in  water 
and  the  value  predicted  by  theory  at  that  temperature  is  an 
indication  of  the  accuracy  of  the  measurement  (approximately 
0.95  percent).  The  density  of  the  sand  was  measured  by 
weighing  1000  cm3 of  water- saturated  sand  and  calculating  its 
mass  The  measured  values  are  shown  in  Table  I . 


111 •  SOUND  PROPAGATION  IN  A  MODELED  SHALLOW  WATER 

ENVIRONMENT 

A.  INTRODUCTION 

The  purpose  of  this  experiment  was  to  investigate  exper¬ 
imentally  the  behaviour  of  sound  propagation  in  a  laboratory 
facility  consisting  of  a  layer  of  water  overlaying  a  layer 
of  sand  and  to  compare  the  results  with  the  normal  mode 
theory . 


B.  EQUIPMENT 

1.  The  tank  used  for  experimental  work  is  made  of  alumi¬ 
nium  plate  one  centimeter  thick.  The  interior  dimensions  of 
the  tank  are  : 

Length  =  9.95  meters 

Width  =  0.94  meters 

Depth  =  0.24  meters 

The  tank  has  a  carriage  which  can  be  used  to  adjust 
horizontal  and  vertical  positions  of  receiver.  The  rails  on 
which  the  carriage  moves  were  leveled  to  ensure  that 
receiver  remained  at  a  constant  depth  as  it  moved  along  the 
tank.  The  leveling  arrangement  was  not  precise  and  it  was 
found  that  the  total  vertical  variation  of  the  receiver 
was  within  +  1  mm.  All  the  joints  of  the  tank  were  sealed 
with  silicon  sealer  and  the  interior  of  the  tank  was  coated 
with  three  layers  of  primer  paint  (zinc  chromate)  to  avoid 
corrosion.  Wooden  wedges  are  fixed  on  the  inner  sides  of  the 
tank  to  minimize  the  reflection  of  sound.  The  tank  was 
filled  with  water  saturated  sand  to  form  an  11  cm  thick 
bottom  layer  and  a  thin  layer  of  water  (2  to  3  cm)  was  added 
above  the  sand.  The  sand  was  degassed  by  circulating  water 
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Figure  3.6 


Pressure  variation  as  function  of  receiver 
depth: f  =  35  kHz  n  =  1 . 


TABLE  VI 


Sound  pressure  as  function  of  depth: 
f  =  50  kHz ,  n  =  1,  L  =  1.5m 


No  . 

Depth ( cm) 

Volts 

Normal 

1 

0.0 

0.04 

0.06 

2 

0.5 

0.21 

0.31 

3 

1.0 

0.41 

0 . 61 

4 

1.5 

0.58 

0.87 

5 

2.0 

0.67 

1.0 

6 

2 . 5 

0.65 

0.97 

7 

3.0 

0.58 

0.87 

8 

3 . 5 

0.435 

0.65 

9 

4.0 

0.34 

0.51 

10 

4.5 

0.245 

0.37 

11 

5.0 

0.15 

0.22 

12 

5.5 

0.1 

0.15 

13 

6.0 

0.063 

0.09 

14 

6.5 

0.05 

0.07 

15 

7.0 

0.05 

0.07 

16 

7.5 

0.047 

0.07 

17 

8.0 

0.027 

0.04 

18 

8.5 

0.02 

0 . 03 

19 

9.0 

0.04 

0 . 06 
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Figure  3.7 


Pressure  variation  as  function  of 
depth: f  =  50  kHz  n  =  1. 


TABLE 

VII 

Sound 

pressure  as 

function 

of  depth: 

f 

-  70  kHz,  n 

=  1,  L  = 

2.0  m 

o . 

Depth ( cm ) 

Volts 

Normalized 

1 

0.0 

0.095 

0.05 

2 

0.5 

0.53 

0.28 

3 

1.0 

1.15 

0.6 

4 

1.5 

1.7 

0.89 

5 

2.0 

1.9 

1.0 

6 

2 . 5 

1 . 75 

0.92 

7 

3 . 0 

1.2 

0 . 63 

8 

3 . 5 

0.7 

0.37 

9 

4.0 

0.24 

0.13 

0 

4.5 

0.205 

0.11 

1 

5 . 0 

0.2 

0.1 

2 

5 . 5 

0.195 

0.1 

3 

6.0 

0.13 

0.07 

4 

6.5 

0.043 

0.02 

5 

7.0 

0.06 

0.03 

6 

7 . 5 

0.065 

0.03 

7 

8.0 

0. 1 

0.05 

8 

8.5 

0.105 

0.055 

9 

9 . 0 

0.073 

0.04 

2 .  Sound  Pressure  as  a  Function  of  Range 

To  measure  the  variation  of  peak  pressure  amplitude 
with  range,  the  transmitter  was  kept  at  constant  depth  and 
the  receiver's  depth  was  adjusted  slightly  at  each  range  to 
measure  peak  amplitude  at  that  range.  Precautions  had  to  be 
taken  to  ensure  that  the  water  surface  is  disturbed  as 
little  as  possible  because  variation  in  water  depth  had  a 
significant  effect  on  the  pressure  amplitude.  Measurements 
were  made  for  three  frequencies  between  the  cut-offs  of 
the  first  two  modes.  The  data  are  tabulated  in  Tables  VIII, 
IX,  and  X  and  is  plotted  in  Figures  3.9,  3.10,  and  3.11 

respectively.  Assuming  a  simple  model  for  which  the  loss  in 
the  bottom  effects  that  portion  of  the  energy  density  of  the 
normal  mode  found  in  the  bottom,  which  then  "robs"  the 
energy  density  in  the  water  layer,  then  the  effective 
absorption  coefficient  can  be  estimated  as  (Ref.  personal 
communication  with  Prof.  Alan  B.  Coppens) 

a  ^  -L  FAfi  d  B/m 

4  A  7 


Where  Ab/A  is  fraction  of  the  energy  density  of  the  normal 
mode  in  the  bottom  and  1/4F  (F  is  the  frequency  in  kHz) 
yields  the  expected  attenuation  to  be  found  in  the  sand  used 
for  the  bottom  [Ref.  4],  Comparison  of  the  measured  data 
with  the  simplestic  model  is  very  good,  suggesting  that 
indeed  the  sound  field  is  subjected  to  minimal  effects  from 
the  side  walls. 
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TABLE  VIII 


Sound  pressure  as  function  of  range: 
H  =  2  cm  f  =  50  kHz  n  =  1 


Di stance ( m ) 

Volts 

Normaliz 

1.00 

0.45 

1.0 

1.25 

0.34 

0.756 

1 . 50 

0.255 

0.567 

1.75 

0.195 

0.433 

2 . 00 

0.2 

0.444 

2.25 

0.168 

0.373 

2 . 50 

0.108 

0.24 

2 . 75 

0.155 

0.344 

3.00 

0.105 

0.233 

3.25 

0.078 

0.173 

3 . 50 

0.022 

0.049 

3 . 75 

0.047 

0.104 

4.00 

0.059 

0.131 

4.25 

0.057 

0.127 

4.50 

0.043 

0.095 

4.75 

0.042 

0.093 

5.00 

0.031 

0.069 

Sound  pressure  variation  as  function  of  range:  f  =  50  kHz 


TABLE  IX 


Sound  pressure  as  function  of  range: 
H  =  2  cm  f  =  70  kHz  n  =  1 


Distance (m) 

Volts 

Normaliz 

1.00 

2.08 

1.0 

1.25 

1 . 6 

0.769 

1 . 50 

1 . 7 

0.817 

1.75 

1 . 0 

0.481 

2.00 

0.93 

0.447 

2.25 

1 .  1 

0.529 

2 . 50 

0.85 

0.409 

2 . 75 

0.68 

0.327 

3.00 

0.87 

0.418 

3.25 

0.62 

0.298 

3.50 

0.17 

0.082 

3 . 75 

0.61 

0.293 

4.00 

0.58 

0.279 

4.25 

0.36 

0.173 

4.50 

0.315 

0.151 

4.75 

0.3 

0.144 

5 . 00 

0.29 

0.139 

TABLE  X 


Sound  pressure  as  function  of  range: 
H  =  2  cm  f  =  100  kHz  n  =  1 


Distance (m) 

Volts 

Norraaliz 

1.00 

1.95 

1.0 

1.25 

1 . 4 

0.718 

1 . 50 

1 . 4 

0.718 

1.75 

1.05 

0.538 

2 . 00 

1 . 0 

0.513 

2 . 25 

0.88 

0.451 

2 . 50 

0.83 

0.426 

2 . 75 

0.85 

0.436 

3.00 

0.74 

0.379 

3.25 

0.63 

0.323 

3 . 50 

0.7 

0.359 

3 . 75 

0.28 

0.144 

4.00 

0.65 

0.333 

4.25 

0.5 

0.256 

4.50 

0.41 

0.21 

4.75 

0.49 

0.251 

3  .  Sound.  Pressure  Variation  Across  the  Width  of  the 
Tank 

Sound  pressure  variations  were  measured  across  the 
width  of  the  tank  at  a  one  meter  distance  from  the  source 
for  three  different  frequencies.  The  measured  data  shown  in 
Table  XI,  XII,  and  XIII  and  plotted  in  Figures  3.12,  3.13, 
and  3.14.  The  results  show  slight  effects  of  interference 
either  from  the  walls  or  from  the  bottom  of  the  sand  laver. 


Transverse  variation  of  sound 


gure  3.14  Transverse 
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F.  SUMMARY  AND  CONCLUSION 

In  the  preceding  sections,  the  experimental  curves  for 
the  pressure  amplitude  as  a  function  of  depth  and  range  for 
three  different  frequencies  for  which  only  the  first  normal 
mode  propagates  have  been  compared  with  normal  mode  theory 
and  a  simple  model  for  attenuation  of  the  normal  mode.  The 
depth  dependent  curves  both  in  the  water  and  in  the  sand 
agrees  satisfactorily  with  the  theory.  The  range  dependence 
obtained  in  the  water  exhibit  excellent  agreement  with  the 
simple  modeling  of  losses  resulting  from  the  relatively 
highly  attenuating  bottom.  Measurements  both  in  range  and 
transverse  dimensions  suggests  that  the  sound  field  is 
subjected  to  minimal  effects  from  the  side  walls.  The 
leveling  arrangement  needs  improvement  because  at  higher 
frequencies  minor  alterations  in  water  depth  changes  the 
pressure  amplitude  considerably. 
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